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STABLE  VISCOSITIES  AND  SHOCK  PROFILES 
FOR  SYSTEMS  OF  CONSERVATION  LAWS 


snail  amplitude  shock  wave  solutions  of  a  nonlinear  hyperbolic  system  of 

conservation  laws  u  +  f(u)  “  0  are  realized  as  limits  of  traveling  wave 

t  x 

solutions  of  a  dissipative  system  ufc  +  f(u)x  “  v^DiUx^x  +  *** 

+  v"(d  u*”*)  .  The  set  of  such  "admissible"  viscosities  includes  those  for 
n  x 

which  the  dissipative  system  satisfies  a  linearized  stability  condition 
previously  investigated  in  the  case  n  =  1  by  A.  Majda  and  this  author. 

When  n  -  1,  we  also  establish  admissibility  criteria  for  singular 
viscosity  matrices  D  ^ (u) ,  and  apply  our  results  to  the  compressible  Navier- 
Stokes  equations  with  viscosity  and  heat  conduction,  determining  minimal 
conditions  on  the  equation  of  state  which  ensure  the  existence  of  the  "shock 
layer"  for  weak  shocks. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  aquations  of  mathematical  physics  take  the  form  of  nonlinear 
hyperbolic  systems  of  conservation  lavs.  With  small  dissipative  effects 
neglected,  typically  smooth  solutions  must  develop  discontinuities  (shocks)  in 
finite  time.  Re-incorporating  dissipation  helps  select  those  discontinuities 
which  are  physically  relevant.  For  this  purpose,  many  different  sorts  of 
dissipation  will  do;  in  particular,  the  physical  viscosity  is  typically 
degenerate  and  not  convenient. 

In  this  paper provide  an  understanding  of  what  high  order  viscosity 


terms  smooth  the  physical  discontinuities.  A  natural  class  of  "admissible* 

viscosity  terms  is  determined  based  on  a  simple  linearized  stability 

— - 

criterion.  In  addition,  jmw  determine  a  class  of  degenerate  second  order 
viscosity  terms  of  physical  type  which  are  admissible.  These  results  are 
applied  to  the  equations  of  compressible  fluid  dynamics,  to  determine  what 


conditions  ensure  the  existence  of  the  "shock  layer  with  viscosity  and  heat 
conduction.  This  should  be  of  interest  to  others  interested  in  general 
equations  of  state  for  compressible  fluids,  such  as  those  investigating  phase 
transitions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
stannary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


STABLE  VISCOSITIES  AMD  SHOCK  PROFILES 
FOR  8 Y STEMS  OF  CON8ERVATZOM  LAWS 

Robert  L.  Peso 


<1.  Introduction 

Consider  •  hyperbolic  systea  of  ■  conservation  laws  in  one  space  distension, 

(1.1)  ut  ♦  f(u)x  -  0  ,  u  S  ■" 

The  problem  we  consider  is  to  determine  those  matrix  n-tuples  (D1,...,Dn>  with  the 
following  property!  A  simple  discontinuous  solution  of  (1.1)  in  the  fora 

u  x  <  at 

(1.2)  u(x,t)  - 

“r  *  *  •* 

is  the  Halt  of  smooth  traveling  wave  solutions  uV  *  U(~~~~~)  of  an  "approximating"  systea 
of  higher  order, 

(1.3)  u_  ♦  f(u)  -  v(D,u  )  ♦  ...  ♦  vn(D  v<B>) 

t  x  1  x  x  n  x 

as  u  +  0  if  (and  only  if)  the  solution  (1.2)  satisfies  a  suitable  entropy  condition  (sea 
Section  3).  Such  an  n-tuple  is  called  admissible.  A  solution  (1.2)  is  called  a  shock  wave 
if  the  entropy  condition  is  satisfied.  (1.2)  is  a  weak  solution  of  (1.1)  precisely  when 
the  Ranklne-Sugonlot  Jump  conditions  are  satisfied! 

(1.4)  f(uR)  -  flu^)  -  slu^  -  uL)  -  0. 

The  traveling  wave  solution  p('  ••  )  of  (1.3)  is  called  a  shock  profile,  with 
C  ■  (x-at)/v,  o(E)  is  required  to  satisfy  ths  systsa  of  ODEs 

.n 

(1.3)  f(0)  -  f(uT)  -  s(0  -  «.)  -  D,  “  ♦  ...  ♦  2-2 

L  L  >  at  n  „.n 
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together  with  boundary  conditions 


(1.6) 


o(C)  ♦ 

“l 

UA 

as 

as 

c  ♦  — 

c  ♦  - 

♦  o 

as  C 

♦  tm 

,  j  ■  1#.. 

dCJ 

Dj  may  be  a  smooth  function  Dj  (u.u^, . . .u*n" 1  * )  for  j  -  1,...,  n. 

A  system  (1.3)  may  naturally  be  associated  with  (1.1)  in  several  ways.  Physically, 
(1.1)  often  arises  as  a  model  for  a  system  with  small,  high  order  viscosity  and/or 
dispersion  terms.  Prototype  examples  are  the  compressible  Navler-Stokee  equations  in  one 


space  dimension,  and  the  KdV-Burgers  aquation  u 


(u  /2)  +  v  u  +  u  u 

x  xx  xxx 


Secondly, 


weak  solutions  to  the  Cauchy  problem  for  (1.1)  are  not  unique,  and  one  hopes  to  identify 
unique  solutions  matheswtically  as  limits  of  solutions  of  some  regularised  equation.  High 
order  terms  may  be  associated  with  (1.1)  in  a  third  wayt  given  a  finite  difference 
approximation  to  (1.1),  it  often  approximates  to  better  accuracy  solutions  of  an  equation 
with  additional  dissipative  and  dispersive  terms  (S] . 

We  assume  that  the  system  (1.1)  is  strictly  hyperbolic,  so  that  if  A(u)  “  3f/3u  is 
the  Jacobian  matrix,  A(v)  has  m  distinct  real  eigenvalues,  ordered  X  (u)  <  ...  <  X  (u) 

1  H 

with  corresponding  right  and  left  eigenvectors  (u)  and  L^tu) ,  k  -  1,...,m,  satisfying 

0  (A* 


(A  -  X^)r^(u) 

£  •  it  ■  6 

*k  rj  °kj 


V  *k,u)  " 0 


An  eigenvalue  X^(u)  is  called  genuinely  nonlinear  (reap,  linearly  degenerate)  if 
VX^*r j(u)  does  not  vanish  (resp.  vanishes  Identically). 

The  problem  above  originated  with  Gelfand  [2] ,  who  suggested  that  the  entropy 
condition  singles  out  those  simple  discontinuities  (1.2)  which  are  limits  of  traveling  wave 
solutions  of  parabolic  systems  associated  with  (1.1)  (the  case  n  «  1  here).  In  more 
concrete  form,  the  investigation  of  the  "shock  layer"  in  gas  dynamics  dates  back  much 
further  (see  [3]).  Host  previous  work  on  the  problem  has  been  for  the  case  n  •  1»  however, 
Shapiro  [9]  and  Smoller  and  Shapiro  [10]  have  obtained  some  results  in  the  case  n  «  2, 
assuming  genuine  nonlinearity. 
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Tha  present  work  is  based  on  the  analysis  (for  the  csss  n  *  1)  of  Hajda  and  Pa 90  [6J , 
who  describe  a  natural  algebraic  condition  on  the  viscosity  matrix  D  «  D1  called  strict 
stability,  and  show  that  any  strictly  stable  matrix  is  admissible  for  all  weak  shocks 
(lu^-Ugl  small).  They  also  obtain  conditions  which  characterise  (up  to  a  degenerate 
class)  those  matrices  admissible  for  weak  k-ahocks  (those  associated  with  a  particular 
Xk,  see  J  3). 

ibis  paper  extends  the  analysis  of  [6]  in  two  directions.  First,  admissibility 
criteria  and  a  notion  of  strict  stability  are  developed  for  n-tuples  (d^,  ...,  Dn)  for  any 
n  (sections  2  and  3).  Second,  admissibility  criteria  are  established  in  the  case  n  ■  1 
for  singular  viscosity  matrices  D(u),  typical  in  physical  systems  (section  4).  Indeed,  in 
the  last  section  we  apply  our  results  to  the  coatpresalble  Havier-Stokes  equations  of  gas 
dynamics,  determining  minimal  conditions  the  equation  of  state  must  satisfy  to  ensure  the 
existence  of  the  shock  layer  for  weak  shocks,  and  to  ensure  that  the  stability  condition 
holds. 

Following  (6],  the  notion  of  stability  for  an  n- tuple  (D1,...,0n)  may  be  motivated  as 
follows<  Linearize  (1.3)  at  a  constant  State  uQ,  obtaining 

(2.1)  u  +  A(u  )u  «  vD,u  ♦  •••  ♦  v*D  u1""1’ 

c  ox  1  xx  n 

If  (1.3)  la  to  be  a  good  approximation  to  (1.1),  a  reasonable  requirement  to  be  imposed  is 

that  the  Cauchy  problem  for  (2.1)  be  strongly  well  posed.  Independent  of  v  as  v  ♦  0  . 

2 

In  b  ,  using  the  Fourier  transform  gives  this  notion  an  algebraic  interpretation! 
Definition.  He  call  the  n-tuple  (D1,...,Dn>  stable  for  (4.1)  if  for  each  T  >  0  there 
exists  C(T)  such  that 

(2.2)  sup  I  exp  t(-i£A(u  )-v£2D.  ♦  ...  ♦  vn(iC)n+,D  )l  <  C(T) 

<Kt<T  0  1  " 

V>0 

£  real 

He  denote  by  Sn ( u0 )  the  set  of  stable  n-tuples,  considered  topologically  as  a  subset  of 
2 

An  n-tuple  in  the  interior  of  Sn(u0)  is  called  strictly  stable  at  u0. 
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Remark .  Introduce  the  matrix  polynoadal 

P(C>  -  -lU(u  )  -  C2D.  ♦  ...  ♦  (1C)B+1D„ 

O  I  It 

Hie  condition  (2.2)  ie  equivalent  to  the  condition 

(2.3)  sup  ■  exp  tf(t)  I  <  C 
t>0 

C  real 

This  section  is  devoted  to  describing  the  structure  of  the  set  of  stable  n-tuples. 
However,  we  point  out  that  a  aajor  objective  of  this  paper  is  to  prove  the  following! 
Theoree  2.1  Suppose  the  n-tuple  (Dj,...,D  )  is  strictly  stable  at  uQ.  If  n  is  even, 
also  assume  Dn  is  nonsingular.  Then  (D,, .. .,Dn)  is  admissible  for  all  shocks  in  some 
fixed  neighborhood  H  of  uQ.  That  is,  if  u^  and  are  in  N  and  satisfy  (1.4), 

then  a  corresponding  shock  profile  solution  of  (1.5),  lying  in  H  and  satisfying  (1.6), 
exists  if  and  only  if  Liu's  strict  entropy  condition  (R)  (see  J3)  is  satisfied  by  the  jump 
(1.2). 

This  theorem  is  a  corollary  of  Theorems  2.3  and  3.1  below,  we  state  here  another 
corollary  of  Theorem  2.3,  giving  a  convenient  sufficient  condition  for  strict  stability! 
Corollary  2.2.  An  n-tuple  (D1,...,0|J)  is  strictly  stable  at  uQ  if  there  is  a  positive 
definite  symmetric  matrix  R  such  that  EA(uQ)  is  symmetric  and 

i)  RDj  is  symmetric  if  j  la  even 

ii)  BDji^-1  is  positive  definite  if  j  is  odd 

If  n  la  even,  we  also  require  that  Dn  has  distinct  eigenvalues. 

The  basic  result  of  this  section  is  Theorem  2.3  below,  which  characterizes  strictly 
stable  n-tuplea.  Theorem  2.4  completes  the  description  of  the  set  of  stable  n-tuples  for 
odd  n.  The  difficulties  encountered  in  trying  to  extend  the  result  to  even  n  are 
analogous  to  those  Involved  when  (1.1)  is  not  strictly  hyperbolic. 

Theorem  2.3  The  n-tuple  (D1,...D|1>  is  strictly  stable  at  uQ  if  and  only  if  the 
following  conditions  holdi 

i)  >  °»  k-  1,...,m. 

il)  If  C  ♦  0,  then  P(5)  has  no  purely  imaginary  eigenvalue. 

ill)  a)  If  n  is  odd,  the  eigenvalues  of  B^"”1  have  positive  real  part. 
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b)  If  n  is  tvtn,  the  eigenvalues  of  Bft  are  real  and  distinct,  and  If 

m  m 

and  denote  corresponding  left  and  right  eigenvectors 
(with  t.  «r.  “  4.  )  then  1  "  B  ln  2  >  0. 

•  j  *  j  X  A*  1  X 

Th worse  2.4  The  set  Sn(u0>  of  stable  n -tuple*  is  the  closure  of  Its  Interior  if  n  is 
odd* 

In  the  rest  of  this  section,  we  prove  2. 2-2.4.  We  begin  by  developing  necessary 
criteria  for  stability.  If  (2.3)  holds,  then  the  eigenvalues  of  P(  £)  oust  have 
nonpositive  reel  pert  for  all  reel  f; .  Using  this  principle,  we  way  establish: 
Proposition  2.5  Aasune  (Dj,...,Dn)  is  stable  at  uQ.  Then 

i)  i)t01rk(uo)  >  0,  k  -  1,...,* 

ii)  For  any  eigenvalue  k^(5)  of  »((),  Be  «^(£)  «  0,  j  -  1, ...,*. 

ill)  a)  If  n  la  odd,  the  eigenvalues  of  Dnin''  have  nonnegative  real  part. 

b)  If  n  is  even,  the  eigenvalues  of  Dn  are  real,  and  if  they  are 

distinct,  then  l.  Or"  i™'2  >  0,  k  «  (,...,■  (notation  as  in  2.3  iii  b) . 
k  hr 

Proof.  11)  is  IssMdiate.  For  convenience,  we  define 
B(9)  *  -P(tan8)(cos9)n/tan8 

-  (cos8)ni  A(u)  ♦  sin8(cos8)n’V  ♦  ...  ♦  (sln6)ni"~\> 
o  i  n 

Tram  ii)  and  continuity,  the  eigenvalues  m^(0)  of  B(6)  satisfy 
(sin8 )Re  u ^(8 )  >  0  ,  -s/2  <  8  <  s/2  ,  j  -  1, ...,«. 

Setting  8  ■  ♦  s/2  we  obtain  lit  a)  and  port  of  iii  b) .  For  i),  8(0)  -  i*  has  distinct 
eigenvalues,  so  for  saall  9  there  exist  seooth  eigenvalues  8)  and  eigenvectors 
1^(8),  with  Ufc(0)  -  iAk(»0).  \10)  -  *^(0),  satisfying 

(B(9>  -  yk(8))j^(8)  -  0  . 
fferentlate,  set  8  »  0  and  dot  with  *^(»0),  obtaining 

VlW  "  wk(0> 

whence  1).  For  ill  b),  a  stellar  procedure  played  at  8  -  s/2  yields 


u;(s/2) 


Proof  of  2.3»  The  necessity  of  tho  condition*  la  anally  established  for  the  most  part,  by 
con*  1  daring  acalar  parturbatlona  of  °vDn'°n-1  **  appropriate.  To  show  that  DR  must 
have  dlatinct  eigenvalue*  when  n  la  even,  we  raaiark  that  a  Jordan  block  for  a  aingle 
multiple  eigenvalue  My  be  perturbed  In  the  (1,2)  and  (2,1)  poaltiona  ao  aa  to  give  rlaa  to 
coaplax  eigenvalue*. 

The  aufficiency  of  the  condition*  la  eatabliahed  a a  for  the  caae  n  ■  1  in  (6) ,  using 

the  Krelas  matrix  theorem,  and  the  fact  that  B(  6)  My  be  amoothly  diagonalized  for  6 

near  0  (and  for  6  near  +  r/2  if  n  la  even). 

Theorem  2.4  follows  directly  from  2.5  and  2.3,  For  if  n  is  odd  and  (D1,...,D||)  is 

stable,  then  it  is  easy  to  check  that  (D,+6l,D,, ...,D  a,D  +  in  1 6i )  is  strictly  stable 

i  2  n- 1  n 

for  any  6  >  0. 

Proof  of  Corollary  2.2.  Observe  that  if  N  la  any  real  symMtrlc  matrix,  and  z  a 

complex  vector,  then  z*Ms  is  real.  Also  if  M  is  positive  definite,  but  not  necessarily 

2 

symmetric,  then  Re(z*Ms)  >  C|z|  .  Now  suppose  (B(6)  -  u  (8))z  ■  0.  Then 

Re  w  (9)  •  z*Ez  »  (cos6)n  I  (tan8)3  Re  z*ED,i3-1z  . 

3  J  odd  3 

So  for  -n/2  <  9  <  x/2.  Re  1^(9)  +  0,  so  ii)  of  2.3  holds.  Also,  for  9  >  0  small. 

Re  11^(9)  >  C  9,  ao  (0)  -  ^D^r^  >  0,  and  i)  holds.  Similarly,  if  n  is  even.  Hi  b) 

holds,  for  then 

Re  II  (9)  z*Ez  -  (sin9)"  \  (cot8)n  3  Re  z*ED  l3  'z 

3  j  odd  3 

>  C(»/2-9)  for  w/2  -  9  >  0  small. 

i  3  Adulssiblllfcy  for  weak  k-ahocks. 

In  this  section  we  characterize,  up  to  a  degenerate  class,  those  n-tuples 
(Df,...,Dn)  which  are  admissible  for  weak  shocks  of  a  particular  family.  As  in  (6) ,  the 
center  Mnlfold  theorem  is  used  to  find  a  trajectory  connecting  critical  points  in  an 
appropriate  system  of  OOEs. 
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We  begin  by  defining  Liu's  strict  sntropy  condition.  First  oonnldnr  tbs  structure  of 


the  Hugonlot  set  of  eolations  of  the  Sankine-Hugonlot  conditions  (1.4).  Fining  u^,  the 
locel  structure  of  this  set  is  well  known  (1].  in  sane  neighborhood  of  uL«  the  poesible 
eolutions  u^  lie  on  n  curves,  uR  ■  uNp) ,  k  “  1,...,n,  passing  through  u^  with 
corresponding  shock  speeds  s  ■  sk(p),  k  •  1,...,n,  satisfying 
uk(0)  -  u^  sk(0)  -  X,^) 

(3.1)  g* (0)  -  r^)  g (0)  -VVVV 

P  "  t^u^)  *  («k(P)~ttL) 

“k 

Liu's  strict  entropy  condition  for  a  k-wave  (1.2)  with  uR  -  u  (pR)  Is  that 
s(l)  ek(p)  >  s  «  sk(pR)  for  p  between  0  and  pR. 

If  k^la)  is  genuinely  nonlinear  and  lu^-u^l  snail,  this  condition  Is  equivalent  to 
Lax's  shock  Inequalities  (see  [6]). 

Theoren  3.1  Fix  uq  8  f  and  k,  1  <  k  <  n.  Assises  X^u)  is  not  linearly  degenerate 
In  any  neighborhood  of  uQ.  Assuan  that  the  n-tuple  (D1,-..,D||)  satisfies  the  following 
nondegeneracy  conditions  at  uQ: 

1)  Dn  Is  nonsingular 

U)  *k°1rk  +  0 

2  ni  1 

111)  -1E(*-X  )(u  )-£  0,  +  ...  +  (1C)  D  is  nonsingular  for  all  real  (to. 

x  o  i  n 

Then  the  following  are  equivalent t 

X»  *kVk(V  >  0 

2)  the  n-tuple  (Dj,...,Dn)  is  locally  adalesible  for  all  k-shocks  In  a  neighborhood 
of  uQ.  That  Is,  there  exists  4  >  0  so  that  If  u^  and  uR  In  Bjfu^)  satisfy  the  jump 
conditions  (1.4)  for  some  s  -  sk(p_),  then  a  shock  profile  lying  In  B.(u  )  exists 

R  0  O 

connecting  u^  to  uR  If  and  only  if  Liu's  strict  entropy  condition  s(S)  is  satisfied. 

theorem  2.1  is  an  lenadlata  corollary  of  this  theoran,  using  2.3.  We  proceed  to  the 
proof  of  3.1.  Our  first  step  is  to  rewrite  (1.5)  as  an  equivalent  first  order  autonomous 
system  of  OOBs.  Introduce  variables  w^  "  uf-if  for  J  •  0,l,...,n-1  and  Introduce  the 
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parameters  v  ■  and  •  as  additional  variable*.  (1.5)  la  now  written,  in  block  fora. 


The  existence  of  a  shock  profile  satisfying  (1.6)  corresponds  to  the  existence  of  a 
trajectory  of  the  system  (3.2)  connecting  the  critical  point  (u^,  0...0,  uL,  s)  to  the 
critical  point  (uR,  0...0,  u^,  s) .  Our  analysis  is  based  on  the  description  of  the  center 
manifold  of  (3.2)  at  the  critical  point  (u0,  0...0,  uc,  X^(uo>). 

Without  loss  of  generality,  assume  uQ  »  0,  X^  ( )  ■  0.  For  convenience,  introduce 
the  colusin  vector  W  »  (w°-v,w*,...,wn"1,v,s),  so  that  W°  -  w°-v.  Then  (3.2)  is  written 

(3.3)  W?  «  T(W) 

For  the  statement  of  the  center  manifold  theorem,  consult  [6J .  To  apply  the  theorem, 

it  suffices  to  describe  two  invariant  subspaces  for  the  linearization  dT  at  the  critical 

point  0i  algebraic  eigenspaces  corresponding  to  eigenvalues  with  zero  and  nonzero  real 

nm  m 

parts,  respectively.  To  calculate  these,  compute,  in  block  form  on  R  x  R  x  R, 

COO 

o 

(3.4)  dT(0)  -  000 

0  0  0 
_  -« 

where  CQ  is  a  block  companion  matrix. 


Sine*  det  Ofl  40  (and  X^tu^)  *  0)  tha  chacactarlatic  aquation  for  dT(0)  say  be  written 

X~’  det(-A+XD>  ...  +  x"d  )  ■  0 
1  n 

Condition  111)  of  3.1  guarantees  that  dT(0)  haa  no  nonsaro  eigenvalue*  with  zaro  real 
part*.  Oondltion  11)  of  3.1  guarantee*  that  tha  zero  eigenvalue  la  semi  simple .  that  la, 
the  algebraic  alganapaca  for  dT(0)  for  tha  eigenvalue  zero  la  equal  to  ker  dT(0).  This 

kernel  la  apanned  by  ■  ♦  2  vectors,  (Ry#  0,  0),  (0,  0,  1),  and  (0,r^,0),  j  -  1 . a. 

Her*  ■  (rfc,  0...0)  «  and  r^  -  r^(0). 

Let  Y  -  ker  dT(0)  and  X  »  range  dT(0).  Then  Y  and  x  are  complementary 
invariant  subspaces  corresponding  to  eigenvalues  with  zero  and  nonzero  real  parts, 
respectively.  Applying  the  center  manifold  theorem,  we  have  (see  (6) ) t 
Proposition  3.2  Assume  that  (0^,...,  Dn)  satisfies  the  nondegeneracy  conditions  1-iii) 
at  u0  -  0  with  X^lu^)  “  0*  Then  there  axists  S  >  0  and  a  Cr  function  (r>2) 

g  t  *  ♦  X  defined  on  Bj(0)  n  y  ■  (y  e  y|  |y|  <  5  }  so  that 

1)  M*  -  {  x  ♦  y  e  I  x  ■  gty)}  is  a  locally  invariant  manifold  for  the  system 

(3.3). 

2)  g(0)  »  0  and  dg(0)  ■  0.  Thus  M*  is  tangent  to  Y  at  0. 

3)  Any  trajectory  of  (3.3:  which  lies  in  Bfi(0)  for  all  £  lies  in  N* .  in 

particular,  critical  points  of  (3.3)  in  Bg{0)  lie  in  M». 

flie  connection  problem  for  (3.2)  is  immediately  reduced  to  one  for  a  scalar  equation 
as  followsi  Define  a  line  in  Y  parametrized  by  y(n)  »  (nR^.u  ,s).  The  curve 
W(n)  *  y(n)  +  g(y(n))  lies  in  M*  while  ly(n)|  <  6,  and  is  Itself  locally  invariant 
for  (3.3),  because  v  and  s  are  constant,  while  g,  mapping  into  X,  is  of  the  form 
g(y(n))  “  (G(n,uL,a),0,0),  G  e  H™*.  Returning  from  w  to  the  (w,v,s>  coordinates  of 

(3.2),  we  find  that  the  system 

o  1 

w  *  w 

e 

(3.5)  1  2 

V  ■  V 

c 

I  .  n-1 

w  ■  D  (f (w°)-f (xr  )-s(w°-u_  )  -  I  O.w^) 

C  n  “  ^  3 
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-  •‘'V*.  -  *  , 


or  w^  »  T(w)  admits  an  invariant  curva 

win.u^.s)  -  <uL+nrk#0...0)  ♦  Gin.u^.s) 

no  long  as  |y(n)|  “  |nr.  I  +  |u.  -  u  I  ♦  Is  -  X.(u  )|  <  6.  It  follows  from  part  3)  of  tha 

K  la  O  K  O 

abova  proposition  that  tha  point  (u^,  0...0)  lias  on  this  Invariant  curve  if 

u_  6  B,(u  )  and  (1.4)  holds. 

R  oo 

The  flow  on  tha  invariant  curve  w(n,u  ,s)  is  now  determined  by  a  scalar  aquation  for 

la 

n(C), 

(3.6)  •  P{n(uL<a) 

where  T  is  determined  from  the  relation 

w  r(o,u_,a)  -  T  (v(n,u.,s)) 
n  l  l 

The  remainder  of  the  proof,  an  analysis  of  the  connection  problem  for  tha  scalar 
equation  (3.6),  is  virtually  identical  with  that  presented  in  [6]  for  the  case  n  “  1,  and 
is  omitted. 


j4.  Admissibility  for  singular  viscosity  matrices 

As  mentioned  in  the  introduction,  viscosity  matrices  in  physical  systems  are  usually 
singular.  In  this  section  we  establish  quite  general  admissibility  criteria  for  weak 
k-shocks  for  such  singular  viscosity  matrices  D(u)  (in  the  case  n-1 ) .  Our  result  will  be 
applied  In  the  next  section  in  a  physical  example,  the  compressible  Navler-Stokes 
equations. 


In  the  case  n-1,  with  D  -  D^(u),  a  shock  profile  0(5)  must  satisfy  the  system 
(4.1)  D(0)0?  -  f<0)  -  f( 

and  boundary  conditions 


V  -  "U“UI. 


0(5)  ♦  as  5  ♦  -**,  0(5)  ♦  uR  as  5  ♦  " 

Theorem  4 . 1  Fix  uq  e  and  k,  1  <  k  <  a.  Assume  X  (u)  is  not  linearly  degenerate 
in  any  neighborhood  of  uQ.  Assume  the  viscosity  matrix  D  •  D^(u)  satisfies  the 
following  conditions! 


1)  D(u)  has  constant  rank  in  a  neighborhood  of  uQ 
ii)  tkDrk(u<))  +  0 
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Sine*  tot  D  +0  (and  X.  (u  )  «  0)  tha  charactartatlc  aquation  for  dT(Q)  may  ba  writtan 

n  K  O 

A**1  dat(-A+XD,+  ...  ♦  n’t  |  ■  0 
i  n 

Condition  ill)  of  3.1  guarantees  that  dT(Q)  haa  no  nonsaro  eigenvalues  with  zero  raal 
parta.  Condition  11)  of  3.1  guarantaaa  that  tha  zaro  eigenvalue  la  aesUalmpla.  that  la, 
tha  algebraic  elgenapace  for  dT(0)  for  tha  eigenvalue  zaro  la  aqual  to  ker  dT(0).  This 
karnal  la  apannad  by  ■  +  2  vactora,  ( Ry ,  0,  0),  (0,  0,  1),  and  (0,r^,0),  j  *  1,...,n. 

Bara  •  (r^,  0...0)  «  *"*  and  r^  »  r^(0). 

Lat  Y  »  kar  dT(O)  and  X  "  ranga  <JT(0).  Then  Y  and  X  ara  complementary 

Invariant  aubapacaa  corresponding  to  aiganvaluaa  with  zaro  and  nonzero  raal  parts, 
raspactivaly.  Applying  tha  canter  nan i fold  theorea,  we  have  (aae  (6)): 

Proposition  3.2  As sum  that  (D^, . . . ,  Dn)  satisfies  tha  nondegenaracy  conditions  1-111) 
at  uQ  »  0  with  A ^ )  -  0.  Than  there  exists  S  >  0  and  a  cr  function  (r>2) 

g  i  Y  ♦  X  defined  on  Bj(°)  n  Y  •  (y  f  y|  lyl  <  6  )  so  that 

D  *•  -  {  x  ♦  y  e  tTm,m< *  |  x  «  g(y)}  la  a  locally  invariant  manifold  for  tha  system 
(3.3). 

2)  g(0)  *  0  and  dg(0)  ■  0.  fliua  M*  is  tangent  to  Y  at  0. 

3)  Any  trajectory  of  (3.3)  which  lias  in  »4(0)  for  all  £  lias  in  M*.  In 
particular,  critical  points  of  (3.3)  in  BjtO)  lie  in  M*. 

Tha  connection  problem  for  (3.2)  la  lamediately  reduced  to  one  for  a  scalar  equation 
as  follows!  Define  a  line  in  Y  parametrized  by  y(n)  “  (nR^.u^.s).  The  curve 
W(h)  “  y(n)  ♦  g(y(n))  lias  in  K*  while  |y(n)l  <  5,  and  is  itsalf  locally  Invariant 
for  (3.3),  because  v  and  s  ara  constant,  whila  g,  mapping  into  X,  is  of  tha  form 
g(y(n))  -  (Gfn.u^.s) ,0,0) ,  G  e  rf* .  Returning  from  w  to  the  (w,v,s)  coordinates  of 

(3.2),  we  find  that  the  system 


(3.5) 


n-1 

I 

3-1  3 
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or  w^  ■  T(w)  admits  an  Invariant  curva 

w^.u^.s)  -  (uI+nr)t,0...0)  +  Gtn.u^.s) 

ao  long  aa  |y(n)l  ■  Inr.  |  *  |u  -  u  I  ♦  |a  -  X.(u  )  I  <  4.  It  follows  from  part  3)  of  tha 

K  L  O  R  O 

above  proposition  that  tha  point  (u^,  0...0)  liaa  an  thla  invariant  curva  if 

u_  e  B,(u  )  and  (1.4)  holda ■ 

R  So 

Tha  flow  on  tha  Invariant  curva  w(n,u^,s)  ia  now  data rain ad  by  a  acalar  aquation  for 

1(C). 

(3.6)  “  Kn.u^.s) 

whare  F  la  datarainad  froa  tha  ralation 

w  F(n,u.s)  -  T  (w( n,u  , a) ) 

V)  L  L 

Tha  remainder  of  the  proof ,  an  analyaia  of  tha  connaction  problaa  for  tha  acalar 
aquation  (3.6),  la  virtually  idantical  with  that  praaantad  in  [6]  for  the  caae  n  ■  1,  and 
la  oalttad. 

t<«  Admissibility  for  alngular  vlacoalty  aatrlcaa 

Ra  mentioned  in  the  Introduction,  viscosity  matrices  in  physical  systems  are  usually 
singular.  In  this  section  we  establish  quite  general  admissibility  criteria  for  weak 
k-shocks  for  such  singular  viscosity  matrices  D(u)  (in  tha  case  n-1) .  Our  result  will  be 
applied  in  the  next  section  in  a  physical  example,  the  compressible  Navier-Stokes 
equations. 

In  the  caae  n  »  1,  with  D  -  Df(u),  a  shock  profile  u(5)  must  satisfy  the  system 
(4.1)  D(0)U5  -  ff(0)  -  ftu^)  -  s(0-uL) 

and  boundary  conditions 

0(5)  ♦  u^  as  5  ♦  0(5)  ♦  u^  as  5  ♦  ■ 

Theorem  4.1  Fix  uq  e  tT  and  k,  1  4  k  <  m.  Assume  X^(u)  is  not  linearly  degenerate 
in  any  neighborhood  of  u0>  Assume  the  viscosity  matrix  0  -  D^(u)  satisfies  the 
following  conditional 

I)  D(u)  has  constant  rank  in  a  neighborhood  of  uQ 

II)  ^Dr^Uj)  +  0 
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ill)  For  all  real  t,  the  matrix  [lttA-X^J-D]  (uq)  la  on*  to  one  on  tlta 
subapaoa  ♦  •  Z^,  where 

*-■{»•  R*|(A-0(u  )v  «  rang*  D(u  )} 

2  It  O  O 

Than  the  following  ara  equivalent i 

1)  l.Dt.  (u  )  >  0 

X  K  O 

2)  D  la  locally  admissible  (or  all  k-ahocks  in  a  neighborhood  of  uQ.  That  la,  there 
exiata  6  >  0  ao  that  if  u^  and  uR  in  Bj(uq)  aatiafy  the  Juaq>  conditlona  (1.4)  for 
aoaw  a  ”  s*(pR),  then  a  ahock  profile  lying  in  Bjtu^)  exiata  connecting  ^  to  t|  if 
and  only  if  Liu'a  strict  entropy  condition  e<*>  la  satisfied. 

The  aaln  steps  in  the  proof  of  thla  theorem  are  the  sane  as  those  for  Theorem  3.1 
above,  or  Theorem  3.1  of  (6).  The  difference  is  that  it  is  a  more  delicate  matter  to 
obtain  an  autonomous  system,  like  (3.2),  to  vhich  the  center  manifold  theorem  nay  be 
applied.  Our  approach  la  to  use  an  algebraic  condition  implied  by  (4.1)  when  D  is 
singular  to  eliminate  some  variables,  than  obtain  an  autonomous  system  for  the  remaining 
variables. 

Aa  before,  introduce  v  ■  u^  and  a  as  variables,  writing  (4. t)  as 

D<o)  U5  -  f(0)  -  f(v)  -  s(U  -  v) 

(4.2)  v^  -  0 


We  motivate  our  elimination  procedure  in  the  caae  that  D  is  constant,  t  linear  and 
a  -  Xr(uq).  Then  (4.21)  la  consistent  only  if  P(A-s)(u-v)  -  0  where  P  is  a 
projection  with  ker  P  •  range  D.  Write  u  -  w+v  for  w  in  ker  P(A-s).  In  order  to 
reduce  (4.21)  to  an  equation  for  w^,  we  should  require  that  0  i  ker  P(A-s)  ♦  range  D 
be  one-to-one.  Mote  that  this  entails  Drk  +  0  and  i^D  +  0,  for 
dim  ker  P(A-s)  m  rank  0. 

Returning  to  the  caae  at  hand,  without  loss  of  generality  we  assume  uq  •  0, 

Xk(uQ)  ”  0.  Let  Zy  -  range  0(0).  Recall  that  D(0)  la  one-to-one  on  z2  from  Hi), 
ao  dim  Z2  -  dim  z,  -  b.  We  may  choose  (inductively  on  dimension)  a  subepace  Zj 
complementary  to  both  Z^  and  Zj,  with  dim  Zj  »  a  -  b.  For  u  sufficiently  asm  11,  we 


i 


ii- 


may  choose  a  smooth  projactlon  P(u)  with  tang*  Z}  and  kernel  ranga(D(u))i  not#  that 
kar  PACO)  -  Zj. 

Si  van  (0,v,s)  in  K2**1 ,  writa  0  ■  u3  -  u2  ♦  v,  whara  u3  ia  in  Zj  and  u2  in 

3  2 

Z2.  Wa  aaak  to  express  u  aa  a  function  of  (u, v, a) ,  uaing  tha  conaiatancy  criterion  for 
(4.21), 

(4.3)  P(U)(f(0>  -  f(v)  -  a(0  -  v))  -  0 

Me  find  it  convenient  to  introduce  isomorphisms 


I,  .  A  ♦  Z2 


*3  >  *“■  *  z3 


2  3  ~ 

and  to  writa  u  •  I  w,  u  «  IjW.  Than  wa  can  apply  to  (4.3),  writing 

h(w,w,v,a)  "  0 

whara  h  s  B*~b  x  *b  x  R*  x  R  ♦  R*”b.  In  block  form,  tha  Jacobian  matrix  of  h  at 


(w,w,v,s)  -0  ia 


dh(0)  -  llj’pACOjlj, 0,0,0] 


ainca  kar  PA(0)  -  z2*  The  firat  coaiponent  ia  an  lsoarerphism  on  rf*,  aince  PACO)  ia 
one-to-one  on  Z}  to  ltaelf.  Thus  tha  iatplicit  function  theorem  applies,  so  that  in  a 
neighborhood  of  0  wa  may  write  w  *  w(w,v,a),  and  indeed  the  total  derivative 
dw(0)  «  [0,0,0].  We  may  axpraaa 

U(w,v,s)  -  I3w(w,v,a)  ♦  IjW  ♦  v 
and  replace  (4.21)  by  the  aquation 


_  30. 
D 


a)  w. 


f(U(w,v,a))  -  f(v)  -  a(U(w,v,a)  -  v) 


By  construction,  the  right  hand  aide  lisa  in  range  (D(U)).  The  matrix 
O  >  Rb  *  range  D  la  an  isoaxirphlsm  at  (w,v,a)  -  0,  so  also  in  a  neighborhood,  aince 
D  has  constant  rank.  Wa  may  find  a  smooth  generalized  inverse 

(D  |^)  *  (w,v,s)  »  if  ♦  Rb 

so  that  (D  )^(o  |~)  “I  in  Rb.  We  have  reduced  (4.2)  to  an  autonomous  system  in 
if  x  if  x  R, 
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(4.4) 


“  (D  "  *<*>  -  a(0(w,v,a) 


v) 


which  we  writ*  W^  -  T(W)  for  th*  variable  w  *  (w,v,s). 

H*  proceed  to  apply  th*  c*nt*r  aanlfold  theorem  to  (4.4)  at  W  -  0.  In  block  fora. 
■Inc*  (3u/3w)(0)  »  Ij  and  UO/3y)(0)  *  I  in  R*,  wa  have 


dT(0) 


(DIj)  A(0)J3 
0 
0 


Condition  111)  of  th*  theoraa  implies  that  dT(0)  haa  no  nonaaro  imaginary  eigenvalues. 
Condition  11)  lapll«a  that  th*  eigenvalue  0  la  aemlslmpl*.  The  kernel  of  dT(0)  la 
■panned  by  a  ♦  2  vector*  (0,0,1).  (0,rj,0),  J  -  1,  ...,a,  and  (R,0,0),  where  I2R  -  rk< 

So,  defining  Invariant  aubapac**  X  »  range  dT(0>,  Y  «  ker  dr(0),  we  apply  th*  center 
aanlfold  theorea  aa  in  |3  to  obtain t 

Proposition  4.2  Assume  that  D(u)  aatlafiea  the  condition*  1-111).  Then  there  exist* 

«  >  0  and  a  Cr  function  (r>2)  g  j  Y  ♦  X  defined  on  Bfi(0)nY  in  tP***'  ao  that 

1)  M*  *  {x+y  s  R  1 (  x  ■  g(y)}  la  a  locally  Invariant  aanlfold  for  (4.4). 

2)  g(0)  -  0  and  dg(0)  -  0.  Thua  K*  la  tangent  to  Y  at  0. 

3)  Any  trajectory  of  (4.4)  which  Ilea  in  B5(0)  for  all  £  lies  In  H». 

Aa  In  (3,  the  connection  problea  for  (4.1)  la  immediately  reduced  to  one  for  a  scalar 
eguatlon<  Define  e  line  In  Y  by  yin)  -  (nR.i^.e).  The  curve  w(t>)  -  y(n)  ♦  g(y(n)>  is 
locally  invariant  for  (4.4),  meaning  it  is  composed  of  solution  curves.  The  curve 
CI(W(n))  la  then  composed  of  solution  curves  of  (4.1).  We  may  write  this  curve 


u(n,uL>s)  -  ♦  nrfc(0)  +  G( n,u^,a) 

where  G  -  l3  w  (nR*g(y(n) ) .i^s)  ♦  lj  g(y(n)).  Rote  that  G(0)  -  0,  dG(0)  -  0. 

If  “r  1b  *nd  aatlafiea  (1.4),  then  u„  ■  U(wR,uL,s)  for  some  wR,  and 

(w^u^.a)  la  a  critical  point  of  T(W),  so  lies  in  M*,  hence  on  w(n).  Thus  for  aoa* 
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»„  ■  «L  *  'i,**  ♦ 

The  flow  on  th«  Invariant  curve  0(n,u. , a)  la  determined  by  a  acalar  aquation 

-  rln.u^.a) 

where 

D(o)  Pln.u^.a)  -  f(n)  -  flu^)  -  sio-u^) 

The  raaalndar  of  tha  proof  ia  identical  to  that  for  nonaingular  D(u),  and  may  bo  found  in 

I«] . 

i5.  Weak  shock  layara  in  compressible  fluid  dynamics 

Haro  wa  use  Theorem  4. 1  to  obtain  vary  weak  condition*  on  tha  aquation  of  atata  In  tha 
coapreaslble  Navier-Stokea  equations  in  ona  space  dimension  which  guarantee  tha  exlstanca 
of  weak  shock  profiles.  He  also  make  a  brief  raniark  concerning  tha  llnaarlsod  stability  of 
these  equations. 

In  Lagranglan  coordinates,  the  equations  are  written  in  conservation  form  as 
Tt  -  Vx  “  ° 

,5',)  *t  *  Px  *  (K>* 

*t  +  (PV)X  “  Vx  +  <f  Vx 

Here  x  la  the  Lagranglan  mass  coordinate,  t  is  time,  t  is  specific  volume,  v  is 

velocity,  p  is  pressure,  8  is  temperature,  S  is  energy  density  per  unit  mama,  and 

It  and  k  are,  respectively,  the  coefficients  of  viscosity  and  heat  conduction. 

2 

£  “  e  +  v  /2,  where  e  is  the  internal  energy  per  unit  mass.  He  assume  that  x  and  6 
determine  the  thermodynamic  state  of  the  material,  and  that  e  and  p  are  given  by 
sufficiently  smooth  equations  of  state,  e  -  e(x,8),  p  -  p(x,8>.  m  and  <  are  positive, 
and  may  also  depend  smoothly  on  t  and  8.  x,  8,  and  p  ars  positive. 

He  assixae  that  the  specific  heat  at  constant  volvsse  is  a  positive  function! 

(5.2)  c(t ,8 )  -  se(T,8)  >  0 

So  8  -  8(x,e)  and  with  u  -  (x,v,s),  (5.1)  may  be  written  in  the  fora 

(5.3)  u  ♦  flu)  -  (D(u)u  ) 

t  X  XX 
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we  shall  presently  that  tha  aquation  u  ♦  f(u)  ■  0  with  diffusion  of  hast  and 
moment  neglected  (w^c«0)  la  strictly  hyperbolic  if  and  only  If 
<5.4)  0  <  -d*/dr|#  eonrt.,Bt  3  «2 

Haro  8  danotas  tha  entropy,  and  8  «  S(t,6).  This  function  la  ralatad  to  a  and  p 
through  tha  Gibba  relations, 

(5.5)  Sd8  -  da  +  pdr 

Our  aain  raault  below  la  that  no  additional  conditions  arc  raqulrad  to  ansura  tha 
existence  of  shock  prof  Ilia  for  weak  shocks.  (The  situation  la  different  for  stronger 
shocks i  saa  (8].) 

Theorem  5.1  Fix  u  ■  (t  ,v  ,&  )  8  It3,  t  >  0»  and  assume  that  (5.2)  and  (5.4)  hold  at 
' . .  '  o  o  o  o  o 

u  .  Then  there  exists  S  >  0  so  that  if  u. ,  u  and  a  satisfy  tha  Rankins  Hugoniot 
°  L  R 

conditions  (1.4)  with  u^,  uR  in  )  -  { u |  |u-uq|  <  4},  than  a  shock  profile  solution 

u(x-st)  of  (5.3)  lying  in  ( u^)  exists  connecting  u^  to  if  and  only  if  Liu's 

strict  entropy  condition  a(K)  is  satisfied. 

The  study  of  the  "shock  layer"  in  compressible  fluid  dynamics  has  a  long  history. 

Host  relevant  here  are  the  results  of  Gllbarg  (3]  and  of  Liu  [4] .  Gilbarg  established  the 
axlstanca  of  shock  profiles  for  shocks  of  any  magnitude,  under  two  additional  conditions  on 
the  equation  of  state: 

(5.6)  0  <  d2p/dr2l  .  _ 

S  constant 

(5.7)  p0(T,0)  >  0 

The  convexity  condition  (5.6)  implies  that  the  eigenvalues  X 1  -  -a  and  1^  »  a  are 
genuinely  nonlinear.  (X2  «  0  is  linearly  degenerate,  see  below.)  In  that  case  tha 
entropy  condition  has  a  simple  form.  Liu  introduced  an  entropy  condition  appropriate  for 
the  nongenulnely  nonlinear  case,  and  showed  that,  with  no  heat  conduction  (k  •  0),  shock 
profiles  exist  for  discontinuities  satisfying  his  sntropy  condition  (see  §3).  This  result 
holds  for  strong  shocks  so  long  as  the  Hugoniot  curves  (see  (3.1))  remain  regular. 

Theorem  5.1  is  proved  by  verifying  conditions  1-11  i)  and  1)  of  Theorem  4.1  for  the 
first  and  third  wave  fields  (k  -  1  and  3).  Discontinuities  associated  with  the  second  wave 
field,  called  contact  discontinuities,  cannot  satisfy  tha  strict  entropy  condition  s(8>. 


The  Jacobian  of  f(u)  takas  tha  forts  (with  p(u)  “  p(r,8(T,s-v  /2) ) ) 


A(t,v,S)  ■ 


Tha  viscosity  matrix  is 


0  -1 
Pt<«)  Pv<ul 

|_ vpt (u)  pfvpv 


0 

pa(o> 
(u)  vp.,(u) 


TD(t,V,« 


-MSt(Tf8)  (U-X)v 


whara  \  “  k/c  >  0.  Tha  information  wa  naad  will  ba  computed  after  performing  a 
convenient  change  of  basis  (slaultaneous  similarity  transformation  of  A  and  D) .  First, 
note  that 

Pt(u)  -  Pt(T,a)  ,  p4(u)  «■  pa<r,e),  p^lu)  -  -vp#(T,a>. 

1  0  0 

0  1  0 

aT(T,8)  v  1 

it  follows  that  «  T  1xD  T  “  diag  (0,u,X)  and 


With 

T  - 


T_1A  T  « 


0 

PT<t,e)  +  pbOt(t,0) 


-1 

0 


0 

Pa(r,a) 


^  0  p+aT(T,8)  0  j 

Tha  eigenvalues  of  A1  are  -a,  0,  and  a.  To  see  this,  compute 
2 

(5.8)  -a  -  p(T,e(T,S))T  -  pT(t,e)  *  PB(T,e)  st(t,S) 
and  observe  that  eT (t ,8)  «  -p  frost  (5.5). 

Condition  1)  of  4.1  is  obviously  satisfied.  We  turn  to  condition  ii). 
(r^»  r2,  r3)  of  right  eigenvectors  of  A^  is 


A  matrix  R  ■ 
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Thus  X,  R  •  R  diag  (-o,0,  a)  and  tha  corresponding  matrix  of  loft  eigenvectors  lo 


(Pt  ♦  P.«t> 
(Pt  ♦  P.V 


Ma  compute 


l1D1r1  -  -  (iw  ♦Xp#(T,o)(ptoT(T,8))/2a 


Proa  (5.5)  and  tha  equality  ST0  “  S0T  on*  aay  verify  tha  tharaodyneaile  identity 

(5.9)  8p8(T,9)  -  p  ♦  #t(t<8) 

Also,  p0(T,8)  •  p^lt.e)  c,  so 

^t0fr1  -  *3Dtr3  ”  <»*a2  *X0cp#2  )  /2a2  >  0 
and  1)  of  4.1  ia  aatisfiad.  For  latar  reference,  we  compute 

(5.10)  *2*-aD1*2  “  “*<Pt<T,e)  *  p#(T,a)aT(T,8))  -  -XpT(r,8) 

It  remains  to  verify  condition  ili)  of  4.1  for  k  •  1  and  3.  Taka  k  •  1. 
z2  *  {  u  «  R3  |  auj  “  )  -  {u  C  ft^MAj+ata  t  rang*  0^ 

To  show  It  (Aj+o)  ♦  la  one-to-one  on  tz2,  it;  aufficaa  to  show  (cf(S.lO))  that 


/  PT(T,8)aa2 
\  a9cp^ 


:i]C 


for  any  complex  a.  A  calculation  aiailar  to  (5.8) ,  using  (5.9),  yields 
(5.11)  pt<t,8)  «  pt(t,8)  -  8cp#2(T,e> 

Therefore,  multiplying  by  dlag  (1,1/a)  from  the  left  and  dleg  (1/8c, 1)  from  the  right, 
it  suffices  to  show  that  iTA^  +  is  noneingular  for  any  t,  where 


But  Dj  la  positive  definite  and  A^  syaawtrlc,  eo  this  is  true.  So  ill)  of  4.1  holds 
for  k  *  1.  For  k  -  3,  replace  -a  by  +a  in  the  argument  above.  This  flniehes  the 
proof  of  Theorem  5.1. 
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Me  conclude  with  a  brief  remark  concerning  the  linearised  stability  condition  (2.3), 

2 

where  P(?)  -  -J(H  •>.  Proa  Proposition  2.5  and  (5.10),  a  necessary  condition  for  (2.3) 
to  hold  is  that  pt(t,9)  <  0.  the  nondegenerate  condition 
(5.12)  0  <  -pT(i,e> 

is  stronger  than  (5.4),  by  (5.f1).  In  fact,  (5.2)  and  (5.12)  are  sufficient  to  inply  the 
linearised  stability  condition,  a  fact  used  by  Matsunura  and  Niahida  in  (7]  to  eatablish 
the  full  nonlinear  stability  of  the  constant  state  for  the  coapresaible  Kavler  stokes 
equations  in  three  space  dimensions,  h  proof  that  the  linearised  stability  condition  holds 
is  similar  to  the  proof  of  Iheorsa  2.3  for  n  >1,  given  in  [6]. 
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